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Recent theoretical and experimental results demonstrate a close connection between the super 
Tonks-Girardeau (sTG) gas and a ID hard sphere Bose (HSB) gas with hard sphere diameter nearly 
equal to the ID scattering length am of the sTG gas, a highly excited gas- like state with nodes only 
at interparticle separations \xjt \ = a; no do ~ «id- It is shown herein that when the coupling constant 
qb in the Lieb-Liniger interaction gBS(xje) is negative and |xi2 1 > inode, the sTG and HSB wave 
functions for N — 2 particles are not merely similar, but identical; the only difference between the 
sTG and HSB wave functions is that the sTG wave function allows a small penetration into the 
region \xi2\ < %node, whereas for a HSB gas with hard sphere diameter Gth.s. = inode, the HSB wave 
function vanishes when \xis\ < fflh.s.- Arguments are given suggesting that the same theorem holds 
also for TV > 2. The sTG and and HSB wave functions for N — 2 are given exactly in terms of a 
parabolic cylinder function, and for N > 2, x no d c is given accurately by a simple parabola. The 
metastability of the sTG phase generated by a sudden change of the coupling constant from large 
positive to large negative values is explained in terms of the very small overlap between the ground 
state of the Tonks-Girardeau gas and collapsed cluster states. 

PACS numbers: 03.75.-b,67.85.-d 



X 



If an ultracold atomic vapor is confined in a de Broglie 
wave guide with transverse trapping so tight and tem- 
perature so low that the transverse vibrational excitation 
quantum is larger than available longitudinal zero point 
and thermal energies, the effective dynamics becomes 
one-dimensional (ID) [l|, 0- 3D Feshbach resonances 
Q allow tuning to the neighborhood of ID confinement- 
induced resonances (CIRs) [HQ, where the ID interac- 
tion is very strong, leading to strong short-range corre- 
lations, breakdown of effective-field theories, and emer- 
gence of highly-correlated iV-body ground states. In the 
bosonic case with zero-range repulsion gs6(xj — xt) with 
coupling constant gs — > +oo, the Tonks-Girardeau (TG) 
gas, the exact A^-body ground state was determined in 
1960 by a Fermi-Bose (FB) mapping to an ideal Fermi 
gas [f| , leading to "fermionization" of many properties of 
this Bose system, as recently confirmed experimentally 
d, 0] • It is now known |^-fl~o{ that the FB mapping is of 
much greater generality; when supplemented by an inver- 
sion and sign change of the coupling constant, it provides 
a mapping between the A^-body energy eigenstates of a 
ID Bose gas with delta-function interactions gsS^Xj — X() 
of any strength [Lieb-Liniger (LL) gas 11] ] and those of 
a spin-aligned Fermi gas. 

Practically all experiments on ultracold Bose gases 
with tight transverse trapping (ID regime) are for the 
case where, in addition to tight transverse trapping, there 
is weak longitudinal trapping by a harmonic oscillator po- 
tential mu! 2 x 2 /2 where x is the longitudinal coordinate. 
In that case the strength of effective ID LL boson-boson 
interactions is characterized by a dimensionless coupling 
constant A = gB/2 3 ^ 2 hujx osc where x osc = y/h/muj is the 



oscillator length. For A — > +oo the exact ground state is 
the trapped TG gas, and for A — s- — oo the same TG state 
is still an exact eigenstate (though now highly excited), 
since both cases map to the ideal Fermi gas [5,, :8]. If A 
is negative and finite, the ground state is very different; 
it is McGuire's collapsed cluster state (bright soliton) 
[l~3j |. However, if A is very large and negative, it was first 
shown theoretically [14H16I ] that the system is metastable 
against collapse to McGuire's cluster state in spite of the 
very strong attractive interactions, and exhibits a strong 
similarity with a trapped ID hard sphere Bose (HSB) gas 
with hard sphere diameter a^.s. > nearly equal to the 
ID scattering length am of the system with A < 0. It is 
even more highly correlated than the TG gas, and hence 
was named [14J the "super Tonks-Girardeau" (sTG) gas. 
By changing A suddenly from large positive to large neg- 
ative values by passing through the CIR, the sTG gas 
was recently created experimentally and shown to have 
properties in agreement with the theoretical predictions 
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Exact solution for N = 2: In the simplest case of two 
bosons with LL interaction gB8(x\ — X2) in a harmonic 
trap potential muj 2 (x 2 + x\)/2, all energy eigenstates and 
eigenvalues can be easily obtained. The wave function for 
the center of mass (cm.) coordinate X = (x\ + x 2 )/2 
is ipc.m. = exp[— X 2 /x 2 sc ] and its energy is -E c .m. = tiuj/2 
assuming that the cm. mode is unexcited. Denote the 
relative wave function by <fi(x) with x — x\ — x<i. The 
excited eigenstates ip u (xi,X2) are most easily found by 
FB mapping 0. The Fermi-Bose (FB) mapping relation 
[5( reduces to 4>b(x) = —sgn(x)c/>F(x)- <j>B satisfies the 
same harmonic oscillator (HO) wave equation as 4>f for 
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i^O, and at x = both satisfy a derivative condition 
which for <fip reads 

<t> F '(0±) = (pg B /fi 2 )<p F (0+) = -(pg B /h 2 )^ F (0-) (1) 

where fi is the effective mass to/ 2. Since cf> F is fermionic, 
it is an odd function, but nevertheless it does not vanish 
at x — 0; instead, it changes sign there, but its derivative 
is continuous. For x ^ 0, <j) F satisfies the free particle 
HO wave equation valid for x ^ supplemented by the 
cusp boundary condition (p}. Two different solutions of 
the HO wave equation having nonzero value and nonzero 
derivative at x = are required, such that the boundary 
conditions of (TIJ are satisfied. These cannot be satisfied 
by the usual Hermite-Gaussian solutions; instead, solu- 
tions which are essentially analytic continuations of these 
to nonintegral quantum number n are required. Since <p F 
is necessarily an odd function of x, it is sufficient to obtain 
a solution y of the differential equation only for x > 0, 
taking <p F (x) = y(q) for x > and <fi F (x) = — y(\q\) for 
x < 0, where x = qx osc . The necessary solution vanish- 
ing and integrable as q — > +oo is a parabolic cylinder 
function D u (q) expressible in terms of confluent hyper- 
geometric series $(a,7; z) which is a sum of two terms, 
one even and the other odd in q 18- 23. 
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This solution vanishes like q u e~ H as q 
verges like \q\~ u ~ 1 e q / 4 as q 
consequence since we require it only for q > 0. The 
relative energy eigenvalues are E IC \ = (y + where 
the nonintegral quantum numbers f are determined by 
the derivative condition (TTJ). This leads to the follow- 



ing transcendental equation for the allowed values [21|: 
T(i — i^)/r(— \v) = —A in terms of the previously 
defined dimensionless coupling constant A. Multiply- 
ing by the signum function so as to obtain the bosonic 
(even) solution <Pb(x) — —sga(x)(f) F (x) one obtains a 
bosonic solution continuous and nonzero at x — but 
with a cusp (derivative sign change) there. The total 
wave function including the cm. factor can be writ- 
ten as ■>pB»(xi,X2) = exp[—(xl + x\) /2xl sc \f v {\q\) where 
f v (\q\) = D v (\q\)e q2 / 4 and q = (x 1 - x 2 )/x osc . The ex- 
ponential growth of e q / 4 as |g| — > oo is cancelled by the 
exponential decrease of D„(|(]r|) like e~ q ^ 4 . 

The first excited state is the one with v — > 1 as | A| — >• oo 
and has only one node at \x\ — X2\— x no d c - For |A| = oo, 
4>f is just the first excited HO state qe~ q I" 1 12 1 with 
a node only at the origin, and for finite negative A, 
^nodc > 0. It is clear from Figs. 5(b,c,d) of [2l[ that 0s is 
almost linear in the region < \x\ < x no d c for all negative 



A, so it is well approximated there by a parabola, which 
can by found explicitly from the leading terms of the hy- 
pergeometric series in Eq. © and the Gaussian prefac- 
tor. By Eq. (|2"j) . fl^ . and the transcendental equation for 
^onefinds0 B «'c{|(j|-f^[l-(f + i)g 2 ) + -- -]} where c 
is a normalization constant. The node is then at X no< ie — 
Oid[1 - (f + j){t S ') 2 + •'•]> which is very close to a iD 
when oid *C x osc , the case when |A| = x osc / \/2air> 1. 

The ground state is very different from the excited 
states, being a collapsed state which is an analog, for the 
trapped system, of McGuire's cluster state [13j. It is an 
even solution also expressible in terms of a Z?„, but one 
whose energy approaches — oo as cjb —> — oo (am — s- 0+); 
see Fig. 5 of [21(. For am — > 0+ it is well approximated 
by iI>bq ~ exp(-|a;i - x 2 \/ai D ) exp[— (x\ + xl)/2xl sc \; in 
fact, this expression satisfies the X12 — > contact condi- 
tion exactly and becomes exact for all (£1,2:2) in both 
limits aiD — > 0+ (total collapse) and aiD — > +00 (ideal 
Bose gas). 

Comparison with trapped hard sphere gas: The recent 



theoretical [14Ml6| and experimental [17] results demon- 
strate a close connection between the sTG gas and a ID 
hard sphere Bose (HSB) gas with hard sphere diameter 
ah.s. nearly equal to am, the ID scattering length of the 
sTG gas. Actually, when all interparticle separations |x| 
are larger than the nodal position x no dc, the sTG and 
HSB wave functions are not merely similar, but iden- 
tical, provided that one sets the hard sphere diameter 
tJh.s. = ^nodc- We prove this first in the simplest case 
N = 2. It is sufficient to restrict ourselves to the re- 
gion x > 2 n odc; the solution for x < — x no d c can then 
be obtained by replacing the argument x by \x\. The 
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parabolic cylinder function D v defined by Eq. @ 
satisfies the correct Schrodinger equation for x > x n odc 
as well as the correct boundary conditions of vanishing at 
M = x n ode and infinity and is normalizable, so it is cer- 
tainly an allowable energy eigenstate for hard spheres of 
diameter x no de- The only remaining question is whether 
it is the ground state for that given hard sphere diam- 
eter. Since Schrodinger's equation for the given region 
and boundary conditions is a well-posed Dirichlet prob- 
lem and the solution given by D u is nodeless in this re- 
gion, it follows from Sturm-Lioville theory that it is the 
hard sphere ground state. At |A| = 00, D v reduces to the 
TG ground state [l2j , which is indeed the ground state for 
hard spheres of diameter a; no do = zero. As |A| decreases, 
this node moves out from the origin, and the solution 
given by D u for x > x n odc continues to be the ground 
state for |A| < 00. The only difference between the sTG 
and HSB wave functions apart from normalization is that 
the sTG wave function allows penetration into the region 
\x\ < ah.s.- This holds for —00 < A < 0, and the penetra- 
tion is small when |A| 3> 1 so that x no de = Qh.s. ~ am- 
Furthermore, the sTG and hard sphere solutions have the 
same energy. We emphasize that when the wave function 
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is extended into the interior region |x| < x no de the resul- 
tant sTG wave function is highly excited relative to the 
collapsed McGuire state due to the effects of the strong 
attraction at x = 0, but when restricted to the region 
\x\ > -Tnodo, it is identical with the ground state of the 
hard sphere Bose gas apart from normalization, and the 
energies of the highly excited sTG state and the hard 
sphere ground state are exactly equal 24(. 

Generalization to N > 2: The ground state in the 
(xi, • • • , xn) space including the interaction points Xj — 
xi is the trapped version of McGuire's cluster state, and 
one expects that for aiD — > 0+ it will be well approxi- 
mated by 
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In fact, the LL contact conditions [ll| arising from the 
delta interactions gs5(xj — xi) are satisfied exactly at 
each collision point Xj = xg. Furthermore, this expres- 
sion becomes exact at all (x\,--- , xjv) in both limits 
aiD — > 0+ (total collapse) and oid — > oo (trapped ideal 
Bose gas), as well as when x osc — > oo. In view of the 
ubiquity of such Bijl-Jastrow pair product wave func- 
tions as exact energy eigenstates of ID Bose and Fermi 
gases 0, E2, l25l - [27j , we expect that the sTG state will 
also be well approximated by such a pair product. The 
obvious generalization of the N=2 sTG wave function to 
arbitrary N is 



i/W(xi, • • • ,x N ) = 
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where U(\qje\) = D„(\q je \)e q ^ /4 and q j£ = (xj-xi)/x osc . 
It also satisfies the contact conditions exactly, and be- 
comes exact at all (x\, • ■ ■ , xn) in both limits aio — > 0+ 
(trapped TG gas, where v = 1 and fv = \<lj — qi\) and 
aiD — > 00 (trapped ideal Bose gas, which maps to the 
trapped FTG gas 2|| 2(|, and where v = and /„ = 1). 

There are accurate numerical solutions for N = 3 
which prove the existence of a gas-like excited state with 
nodes only at a nearest neighbor separation \xji\ = x n ode 
which increases with decreasing |A|, is very close to a\D 
for |A| 3> 1, and goes to zero along with am in the TG 
limit |A| — > 00; see Fig. 4 of 22[. Furthermore, for all 
N > 2 the approximate wave functions Q vanish only 
at \xji\ — a; n odc and become exact both at the collision 
points Xji = and when all \xji\ — > 00, and we expect 
the exact sTG solution to have the same properties. We 
expect that the theorem on the identity of the sTG and 
hard sphere wave functions and energies in the region ex- 
terior to the nodes will also hold for all N > 2 since the 
exact hard sphere ground state for hard sphere diameter 
Xnodc = is the |A| = 00 TG ground state [12[ and the 
same arguments we used for N — 2 should apply as |A| 



is decreased. A rigorous appeal to Sturm-Liouville the- 
ory would require a proof that for the exact, unknown 
sTG eigenfunctions the pair distribution function van- 
ishes only at \xj — xg | = x n qdc; and so far this is known 
only for N — 2 and N = 3 [22j. Therefore we state our 
belief of the identity of the sTG and hard sphere wave 
functions and energies for all N as a conjecture. 

Sudden approximation and metastability: In the ex- 
periment by the Innsbruck group (l7| . the dimensionless 
coupling constant A is changed suddenly from large pos- 
itive to large negative values. For A > 1 the ground 
state is very close to that of the trapped TG gas, the 
ground state for A = +00, whose exact wave function 
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is a normalization constant 



x j /2x osc ) 
12( . Assum- 



is IpTG = 

where C n 

ing that the switch to large negative A is rapid enough 
that the sudden approximation is valid, the wave func- 
tion after the switch will be in a superposition of all 
eigenstates -0a q of the system with negative A: ip\(t) — 
Y Ja ( % l ; ^\' l kTG)^Xae~ lExat/h . The stability of a trapped 
gas with attractive ID interaction was addressed in 
17] by measuring the frequencies of collective breath- 
ing mode. It was found that by rapidly crossing the 
CIR to finite positive am the frequency increases, which 
agrees with the description of the sTG state in terms of 
an equivalent hard rod system, while at larger values of 
am the frequency decreased. These results agree with 
our theorem that the sTG state has the same ground 
state energy as a system of hard spheres, increasing with 
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29f . so that the frequency in the sTG state 
increases with a\o, although for non-zero values of am 
there is a possibility of populating states different from 
the sTG state. We analyze the transition through the 
CIR in terms of the projection of the TG state to the gas- 
like sTG eigenstate and states that contain bound states, 
relying on the fact that the initial state for A 3> 1 is very 
close to the A = +00 TG state. In the above expres- 
sion for ip\(t) there are many-body bound states with 
N = 2, ...,7V particles, with the bound part described 
by McGuire's many-body solution Eq. ([3]) with a char- 
acteristic size of am The projection integral {i^xo^tg) 
decreases exponentially fast with the number of particles 
in the bound state. From this the most probable way 
to transfer particles to some bound state is to do so to 
a bound state with 2 particles. We calculate the projec- 
tion integral (V'AqIV'tg) f° r different numbers of particles 
and show the results in Fig. [T] The transition to many- 
particle bound states is highly suppressed compared to 
the transition to a state where two particles are bound 
while all the other particles are unbound. This can be 
understood my comparing characteristic volumes that a 
state occupies in the phase space. The typical size of a 
McGuire state of ./V particles is a^ D , while the TG gas 
occupies (-\/Nx osc ) N , so that the overlap decreases dra- 
matically as N is increased. From this the most relevant 
transition is to a state with two particles bound. The 



4 




FIG. 1: Square of overlap integral of TG gas-like state with a 
McGuire state containing N c i clustered particles, as a function 
of am/xoac. Inset, the same data on a semi- logarithmic scale. 



overlap integral was evaluated analytically for 2 bound 
particles in systems of N = 2,3,4 particles, and numer- 
ically for bound clusters of 3 and 4 particles, with the 
bound cluster described by the McGuire wave function 
Q while leaving the other particles in the TG state. In 
all cases the overlap is very small when aiD *C x osc , as 
in the experiment [17| . 

These results are expected since ipTG is an exact eigen- 
state not only for A = +oo, but also for A = — oo, since 
both map to the ideal Fermi gas [5, 8]. Hence the system 
is completely stable under a sudden jump from A = +oo 
to A = — oo, and metastable following a jump from A > 1 
to A < -1. 
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ing the proof of the identity of the sTG and hard sphere 
wave functions outside the node at \x\ = x no d e - Research 
of M.D.G. is supported by the U.S. Army Research Lab- 
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